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ANTON LEYKIN 

Abstract. Wc use our recent implementation of a certified homo- 
topy tracking algoritlim to search, for start systems that minimize 
the average complexity of finding all roots of a regular system of 
polynomial equations. While finding optimal start systems is a 
hard problem, our experiments show that it is possible to find 
start systems that deliver better average complexity than the ones 
that are commonly used in the existing homotopy continuation 
software. 



1. Introduction 

This section introduces homotopy methods for solving polynomial 
systems and sets up the stage for the experiments that are in the core 
of this article. The main contribution of this work concerns optimal 
start systems, which do not enter this stage until Subsection 11.21 

The problem of finding a good algorithm to solve a system of poly- 
nomial equations has been on the minds of many mathematicians for 
many centuries. Here we address this problem in its most basic form: 

• the coefficients of the polynomials are complex numbers; 

• the system is square (neither under- nor overdetermined) ; 

• the system has finitely many complex solutions and this number 
is equal to the the number of solutions of a generic system 
of equations with the same degrees (the Bezout bound) and, 
therefore, all solutions are regular. 

In principle, the stated problem can be solved symbolically. Let field 
C be an extension of Q containing all coefficients. Then, for instance, 
via elimination techniques employing Grobner bases one can find a 
finite extension of C that contains all coordinates of all solutions. In 
practice, for systems with a large number of solutions this method 
could be impractical: the complexity of expressions representing the 
resulting extension and solutions is high, not to mention the phenomena 
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of intermediate expression swell. Besides, even if these expressions are 
obtained, getting numerical approximations to the solutions could be 
a nontrivial task. 

Approximate solutions would arguably constitute the most widely 
applicable output of any polynomial system solving algorithm. A basic 
question to ask at this point is: What does it mean for a solution to 
be approximate? Here are two answers: 

• fix e > 0, then z G is an approximation of an exact solution 

C e if ll^-CII < 

• fix a refinement algorithm, i.e.: a function i? : C" — ?■ C", then z 
is an approximate zero associated to ( if the sequence {R''{z)} 
converges to (. 

The first way, while frequently preferred by practitioners, is clearly de- 
pendent on the threshold e. On top of that, one relaxation is commonly 
made in applications: the distance — Cll is replaced by its estimate. 

The second approach is intrinsic once the refinement procedure R 
is fixed. Set R to be Newton's iterator to get the notion of Smale's 
approximate zero used as a basis for a-theory. (These concepts are 
defined in Section [2J) It is this meaning of an approximate solution 
that we adopt for the purpose of this article. 

For the purpose of finding approximate solutions let us introduce the 
homotopy continuation technique: given a target system f construct a 
start system g and track the curve 

t ^ ht, t E [0, T]; such that hr = f and = g; 

in the space of polynomial systems. An atomic task of a homotopy 
continuation algorithm is to take a start solution, i.e., a zero zq G 
g~^{0), and track the homotopy path starting at zq and finishing at a 
target solution, a zero of /. For details on how g and t ^ ht are set 
up so that the last sentence makes sense one may consult the recent 
book [16]; we construct several particular homotopies in Section [2l 

In the last decade several software systems (e.g.: [T], [5], [lU], [T7] ) 
based on homotopy continuation powered by numerical predictor-corrector 
methods have been developed. Problems with millions of solutions 
could be attacked due to speed achieved through path-tracking heuris- 
tics. The target solutions in some cases can be certified rigorously, how- 
ever, usually only data characterizing the quality of a solution (such as 
estimates for condition number, error, residual value, etc.) is produced. 
The latter is aimed at providing the software user with confidence in 
the obtained results, which is arguably what the most users need. 
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1.1. Certified numerical homotopy tracking. Apart from the cer- 
tification of the end points of the homotopy paths, we would also like 
to certify that all points produced by the tracker working on one homo- 
topy path are approximate zeros associated to the points on this path. 
Note that a heuristic tracker could lead to the picture in Figure [1) at 
every point where a prediction is evaluated (depicted by arrows) the 
system is well-conditioned, approximate solutions to all solutions of 
the target system are discovered, however the continuation paths are 
swapped by the tracker. 




Figure 1. "Path-crossing" scenario: the tracker jumps 
from path 1 to path 2 and vice-versa. 

Ensuring this does not happen is of ultimate importance to the al- 
gorithms that produce discrete output based on the assumption that 
no path-jumping occurs. For instance, discovering monodromy via ho- 
motopy tracking is a part of algorithms for the genus of an irreducible 
curve within an algebraic set [2], Galois groups of Schubert problems 
[TT] . and numerical irreducible decomposition. All these are applica- 
tions to problems in pure mathematics and certification is highly desir- 
able in order to give the obtained results the status of a mathematical 
proof. 

Creation of a certified homotopy tracking procedure was the goal of 
our joint work with Carlos Beltran in the resulting algorithm has 
been implemented in the NumericalAlgebraicGeometry package [10] of 
Macaulay2 [8], referred to as NAG4M2. 
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1.2. Start systems with optimal average complexity. We used 
the implementation mentioned above to estimate the average complex- 
ity of computation with the linear homotopy using various initial pairs 
(start system, start solution). The results concerning the average com- 
plexity of finding one solution of a system obtained in |1] are outlined 
in subsection 14.11 

The main question considered in this paper and not addressed in |1] 
is: How to find an optimal start system, i.e., one that delivers the 
minimal average complexity of computing all solutions^Q 

This is a hard question, which has not been answered even in the 
one-equation case. The work of Shub and Smale [13] suggests that an 
optimal start system for the linear homotopy has to have the smallest 
condition number. They define the concept of the condition number for 
a polynomial equation and prove that there is a family of polynomials 
— elliptic Fekete polynomials — for which the condition number grows 
polynomially with the degree. 

Subsection 14.11 revisits the experiments of |1] with a view towards a 
search of an optimal initial pair (for finding one solution of a target 
system). In subsection 14.21 the results of several new experiments are 
presented in detail show that it is possible to find start systems that give 
better average complexity (of finding all roots) than the total-degree 
homotopy systems that are a popular choice in current software. 

While the discovery and experiments are done with the certified 
homotopy tracking algorithm, heuristic computation can benefit from 
finding either optimal or near-optimal start systems as well. This is 
backed up by the data obtained by running an implementation of a 
heuristic algorithm in several experiments. 

The Macaulay2 [5] and Mathematica [IS] scripts used in this paper 



are published at http://people.math.gatech.edu/~aleykin3/0SS/ 
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In the context of structured polynomial systems, the word optimal is used in a 
different sense: it refers to a homotopy that minimizes the number of continuation 
paths tracked in the computation. For example, polyhedral homotopies imple- 
mented in jl71 [3] are optimal for solving sparse polynomial systems in that sense. 
The exact meaning of optimality in the context of this paper shall become clear in 
Section m 
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2. Preliminaries 

To simplify the discussion in this paper we shall always consider 
the homogeneous problem of polynomial system solving. For a positive 
integer / > 1 , let be the vector space of all homogeneous polynomials 
of degree I with complex coefficients and unknowns Xq, . . . , Xn. For a 
list of degrees (d) = {di, . . . , dn) let 

'H(d) = T^dx X ■ ■ ■ X T-Ld^ . 
An element h G 'H(d) is seen as a vector in a space of dimension 

and as a system of n homogeneous equations with n + 1 unknowns. The 
zero set of /i is a projective variety: we consider its zeros as projective 
points C G P(C'"'*"^). Throughout the paper no distinction between a 
point in P(C""'"^) and a representative of the point in C^"*"^ is made: 
when necessary it is implied that a representative has the unit norm. 

2.1. Approximate zeros and the projective Newton's operator. 

It has been pointed out in the introduction that describing the zeros 
of /i G T-i{d) exactly is a hard task: one may ask for points which are 
e-close to some zero, however, this is not an intrinsic concept. What 
we employ here is the concept of an approximate zero of [T5] . 

First we define the projective Newton's method as follows. Let h G 
and z G P(C"+i). Then, 

Nv{h){z) = z ~ {Dh{z) |,x)-'/i(^), 

where Dh{z) is the n x (n + 1) Jacobian matrix of /i at z G P(C"^"^), 
and 

Dh{z) 

is the restriction of the linear operator defined by Dh{z) : C""*"^ — ?■ 
to the orthogonal complement z-^ to the one-dimensional space spanned 
by z. Hence, Np{h){z) is defined as long as {Dh{z) \z^)~^ , a linear 
operator from to z-^ C C^"*"^ of dimension n, is invertible. 
Let dR be the Riemann distance in P(C"^^): 

dR{z,z') = arccos^||J| G [0,7r/2], 

where (■, ■) and || ■ || are the usual Hermitian product and norm in C"~*"^. 
Note that dR is well defined on P(C"+^) x P(C"+^). 
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Definition 1. We say that z G P(C""^"'^) is an approximate zero of 
h G Ti^d) with associated zero ( G P(C""'"-'^) if Np{hy{z) is defined for 
all I > and 

dRiNpihyiz)x)<^^^, i>o. 

2.2. Bombieri-Weyl norm. Given two polynomials v,w E Tii, 

«0 + --- + "n=' 

ao+...+an=l 

we define their Bombieri-Weyl product to be 

-1 



(^5^) / , ( / \ ) ^ao,...,anbaQ,...,any 

ao+ai + ...+an=l 

where ~ is the complex conjugation and 

/ \ /! 



,(ao, • • • ,«n)/ aol---aJ 

is the multinomial coefficient. 

Then, given two elements h = {hi, ... , h^) and h' = {h[, . . . , h'^) of 
T-L^d), we define 

{h, h') = {hi, h[) + --- + {hn, h'^), \\h\\ = ^/ihTh). 

From now on, we will denote by S the unit sphere in 'H(d) for this norm, 
namely 

§ = {/iG?/(d): \\h\\ = l}. 

2.3. The condition number. The (normalized) condition number dX 
{h,z) G Hi^d) X P(C"+i) is defined as follows 



n{h, z) = \\h\\ {Dh{z) l^x) """Diagdl^; 



or 2) = oo if Dh{() \ ^± is not invertible. Here, \\h\\ is the Bombieri- 
Weyl norm of h and the second norm in the product is the operator 
norm . Note that assuming the system and a representative of a pro- 
jective point are normalized, i.e., /i G § and \\z\\ = 1, 

^I{h,z) = \\{Dh{z)\,.)-'Dia.g{^/dl)\\ 

is determined by the operator norm of the inverse of the Jacobian Dh{Q 
restricted to the subspace orthogonal to z multiplied with a diagonal 
matrix that makes formulas look nicer. 
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We also define the condition number of tlie system h as 

(2.1) fi{h) = sup fi{h,z). 

Note tliat witli tliis definition /i is a regular polynomial system if and 
only if fL{h) < oo. 

2.4. A total degree homotopy. For a real number r > (a common 
setting is r = 1), let the start system be 

(2.2) g = {Xf' - r'^'X^' , . . . , X^" - r'^"X^"). 

This system is used in practice by various homotopy continuation soft- 
ware packages due to the simplicity of its evaluation and the simplicity 
of its solutions. 

In NAG4M2 a total-degree start system (with r = 1) is created as 
follows: 



il 
i2 
i3 



R = CC[x,y,z] ; 

f = {x~2+y~2+z~2 , x*y}; — target system 
totalDegreeStartSystem f 



2 2 2 2 
o3 = ({x - z , y - z >, 

{(1, -1, 1), (-1, 1, 1), (1, 1, 1), (-1, -1, 1)}) 

2.5. A conjecture by Shub and Smale. The condition metric on 
the so-called solution variety 

V = {{h, z)eSx P(C"+^) I h{z) = 0} 

is obtained by multiplying the metric inherited from the product S x 
P(C"~'^^) by the condition number fi. The length of a homotopy path 
t — {ht, Ct) then equals 

(2.3) Co(/,^7,Co)= f fi{ht,Ct)\\ht,Ct\\dt. 







Consider the following initial pair: 



(2.4) g={-: , eo 



The following is a slightly modified conjecture of Shub and Smale (the 
pair {g, cq) is a slightly modified version of the pair in the original 
conjecture) : 
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Conjecture 1. [TH Conjecture 2.4] For g picked in the sphere S with 
the uniform probability distribution, the expected length of the homotopy 
path E(Co(/, (?, Co)) is bounded by a polynomial in N . 

Proving polynomial average complexity for the linear homotopy with 
the initial pair above would give a deterministic algorithm settling 
Smale's 17th problem. However, the tightest known bound on this 
complexity obtained so far is A^'-'(i°gi°s^)j see Biirgisser and Cucker [7]. 
Most of the recent developments in this direction are inspired by |12] . 

In NAG4M2 a good initial pair is created as follows: 

11 : R = CC[x,y,z] ; 

12 : f = {x"2+y"2+z"2, x*y}; — target system 

13 : goodlnitialPair f 

o3 = ({1.41421x*z, 1.41421y*z}, {{0, 0, 1}}) 

Note that the start system g of (12. 4p has cq as the only isolated 
solution with the smallest condition number possible, while the rest 
of zeros are described by Xq = and form a set isomorphic to P"~^. 
In particular, system g can't be used to compute all solutions of the 
target system /. 

2.6. Random linear homotopy. A randomized approach has been 
developed by Beltran and Pardo [5] (see also |6]): the basic idea is 
to construct an initial pair [g, zq) in a random fashion. First, pick a 
random start system g in the sphere § with uniform distribution, then 
pick Zq to be a random zero of g. While the first task is straightfor- 
ward, a naive approach to the second would be dependent on solving 
g. Nevertheless, there is a clever way to build a pair with the given 
properties pointed out in [5] that depends only on the ability to pick 
random points in a sphere with the uniform probability distribution. 
The procedure is implemented in NAG4M2 package along the descrip- 
tion given in 

11 : R = CC[x,y,z] ; 

12 : f = {x~2+y~2+z~2 , x*y}; — target system 

13 : randomlnitialPair f 

2 

o3 : ({(- .11245 + .222955*ii)x + (.28313 + . 30427*ii) x*y + ... , 

2 

(.321684 + .0258085*ii)x + (.200917 + . 161032*ii) x*y + ... >, 
{{.30545-.0745817*ii, . 277793+ . 71328*ii , . 550537- . 11004*ii») 
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Thus constructed random linear homotopy is shown to have average 
polynomial complexity in [S] providing a uniform randomized algorithm 
that solves Smale's 17th problem. 



3. Certified linear homotopy 
Here we give a partial summary of the main constructions of |4]. 

3.1. Linear homotopy. Given two systems / and (7 in the unit sphere 
S C 'H(d) we define a linear homotopy as a segment of the geodesic curve 
on S connecting / and g given by its arc-length parametrization: 



(3.1) t ht = g cos{t) 
where 



/-Re((/,^))^ 
y/l-Re{{f,g)y 



sin(t), te[0,T] 



T = arccos Re((/, = distance(/, (?) G [0, 



TT 



The procedure of certified tracking for a linear homotopy is presented 
by Algorithm [TJ 

Algorithm 1. z^, = TrackLinearHomotopy{j\ g, Zq) 

Require: f\g ^ S; zq is an approximate zero of g. 
Ensure: z^ is an approximate zero of f associated to the end of the 
homotopy path starting at the zero of g associated to Zq and defined 
by the homotopy liS. 
1: i ^ 0; Si = 0. 
2: while Si ^ T do 
3: 



Compute 

gi ^ hs 



at s = Si. 

"Pi ^ where 



f-Re{{f,g))g 
v/l-Re((/,^?))2 



cos s 



Dgi{zi 
z* 



\ 



1/ 



Xi,2 



^ \\9i\ 



Dgi{zi 
z* 



gi[Zi 
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5; 



Let ti he any number satisfying 



0.04804448 



0.04804448 



<t,< 



6: if ti > T — Si then 
7: ti^T~Si. 

8: end if 

9: Sj+i ^ Si+ ti; Qi^i ^r- /is,+i- 

10: Perform a step of the projective Newton's method: 



11: i-^i + l. 
12: end while 

13: z^, ^ zt- 

The correctness of the algorithm is shown in [3]. The complexity of 
the algorithm is bounded by the following: 



where d the maximal degree of polynomials in the system and Co is 
defined by (10) . 

The following is an example of usage of one of the main functions of 
NAG4M2 with options that specify the computation to Algorithm [1] 

11 : needsPackage "NumericalAlgebraicGeometry" ; 

12 : f = randomSd {2,2}; — target system (random 2 quadrics) 

13 : (G, solution) = randomlnitialPair f; 

14 : trackCg, f, solution, Predictor=>Certif ied) 

04 = {{.050877+.571108*ii, -. 74469- . 020102*ii , -. 154178- . 304176*ii» 

15 : (first 00) .NumberOf Steps 

05 = 781 

In this example a solution of a random system is found by a linear 
homotopy with a random initial pair in 781 steps. 



First we give outline several several experiments staged in and 
then present the results of an experimental search for optimal start 
systems. 



Zi+i ^ \\Np{gi+i){zi)\\ Np{gi+i){zi). 



projective Newton's method steps | < \71d^^'^Co\ 



4. Experiments 
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4.1. Average complexity comparison. We have obtained experi- 
mental data by running a linear homotopy connecting an initial pair 
{g, zq) to a random system in S C 'H(rf) with di = 2 for i = 1, ... ,n. 
Three kinds of initial pairs were examined: 

• good: The initial pair (12. 4p conjectured to be "good" by Shub 
and Smale; 

• total: The start system g of the total degree homotopy (12.21) 
with zo = (1, !,...,!); 

• random: The random initial pair discussed in Subsection 12.61 

One can see the details of this experiment in [1]: to summarize, the 
conclusion is that the average complexity of good is smaller than that 
of total, which in turn is better than the complexity of random (for 
all n, for which we ran the computation). 

Here we report a stronger experimental negative result: in experi- 
ments with the problems having at most 8 solutions we failed to find an 
initial pair that performed better than the "good" initial pair. (More 
precisely, a one-day long random search for each case with systems hav- 
ing at most 8 solutions fails to get an initial pair that performs better 
than (12.41) on average.) 

This not only provides an experimental evidence to, but also prompts 
the following stronger version of the conjecture by Shub and Smale: 

Conjecture 2. The initial pair [g, cq) as described in ^2.4^ minimizes 
E{Co{f,g, Zo)) over all initial pairs {g, Zq) in the solution variety V. 

The weaker Conjecture [1], indeed, follows from Conjecture [2] in view 
of the existence result [TH Main Theorem]. 

4.2. Optimal start system. Much like in the problem of finding one 
zero, in [13] Shub and Smale use the condition number of a system as 
defined by (12.11) as a guideline for picking a "good" start system. They 
proceed to show that in case of n = 1 (one equation) the condition 
number of the total-degree start system (12.21) grows exponentially as a 
function of the degree d of the equation. In turn, the so-called elliptic 
Fekete polynomials give a family for which fi{Fd) = 0{d). The 
construction of this family goes via Fekete points on the two-sphere, 
the problem of finding which relates to Smale's 7th problem. 

To summarize, minimizing the condition number of a start system is 
a hard problem. The experiments below suggest that the complexity of 
the computation with a given start system, which ultimately is linked to 
the sum of the length of homotopy continuation paths in the condition 
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metric 

E{Co{f,g))= Yl E(Co(/,^7,z)), 

is not determined just by its condition number. 

Our methodology is simple to describe (the choices of constants were 
made in such a way that the runtime of each of the computations 
reported in this article does not exceed one day): 

(1) Among 1000000 random start systems g picked on § with uni- 
form probability distribution we select five with either the small- 
est fi{g) or the best estimated average complexity; 

(2) Approximate the average number of steps the implementation 
of Algorithm [T] takes to compute all solutions by averaging over 
10000 random target systems; 

(3) Compare with the average performance of specially constructed 
start systems (total-degree, elliptic Fekete polynomials). 

For the construction of the total-degree systems with the minimal 
condition number we have created a script in Mathematica that con- 
struct a (symbolic) expression for the condition number ix{g) depending 
on the parameter r in (12. 2 p and performs either symbolic or numerical 
optimization determining the optimal value for r. 

4.2.1. Case n = 1. Consider a single equation of degree four: one can 
look at 

gtotal = — 

as r = 1 in (12. 2p is optimal for this case, and 

gPekete = -^1 (-^1 ~ 2a/2Xq). 

The construction of the latter can be carried out explicitly using the 
general recipe of [13] . Four Fekete points are the vertices of the regular 
tetrahedron inscribed in the sphere S'^ of diameter 1 and the center at 
(0,0, |). Placing one of the vertices in the origin (see Figure [2]) the 
stereographic projections of the others onto the xy-plane lie on a circle 
of radius a/2. 

The results of two searches are reported below. Both use the method- 
ology described in the beginning of the subsection, but differ in the way 
the systems are chosen in step (1). 

• The search presented in Figure |3] is done according to the con- 
dition number ^{g); 

• The search presented in Figure H] uses the routine that estimates 
the average complexity. 
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Figure 2. Construction of an elliptic Fekete polynomial 
from four Fekete points. 



The tables contain the detailed conditioning data as well as the esti- 
mated complexity (average number of steps in the certified homotopy 
continuation algorithm) for the specified start systems. The random 
systems Qi in the tables are ordered according to fi{gi). 





91 


92 


53 


94 


95 


9Fekete 


9total 


1^(9,^2) 

IJ'{9,Z3) 

K9,Z4) 


1.23188 
1.23644 
1.21733 
1.23051 


1.23094 
1.23788 
1.23569 
1.22567 


1.23585 
1.23836 
1.19765 
1.23537 


1.23975 
1.22948 
1.23153 
1.23855 


1.24018 
1.20723 
1.23769 
1.23718 


1.22475 

1.22475 
1.22475 
1.22474 


1.41421 
1.41421 
1.41421 
1.41421 


^9) 


1.23644 


1.23788 


1.23836 


1.23975 


1.24018 


1.22475 


1.41421 


#steps 


1138.32 


1136.93 


1137.81 


1141.76 


1130.81 


1137.51 


1177.79 


Figure 3. Search by condition number: n = 1, d = 4. 




91 


92 


93 


94 


95 


9Fekete 


9total 


K9,zi) 
^(5,^2) 

K9,Z3) 
K9,Z4) 


1.28246 
1.25521 
1.23108 
1.17805 


1.30382 
1.14504 
1.38061 
1.35972 


1.33386 
1.3199 
1.35796 
1.14828 


1.43745 
1.10411 
1.28299 
1.41302 


1.44709 
1.34199 
1.25412 
1.25466 


1.22475 

1.22475 
1.22475 
1.22474 


1.41421 
1.41421 
1.41421 
1.41421 




1.28246 


1.38061 


1.35796 


1.43745 


1.44709 


1.22475 


1.41421 


#steps 


1138.47 


1146.18 


1152.21 


1154.79 


1155.42 


1137.51 


1177.79 



Figure 4. Search by average complexity: n = 1, d = A. 

We can see that both the elliptic Fekete polynomial and discovered 
by a random search perform significantly better than gtotai- A natural 
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question arises that, to our knowledge, has not been answered yet: Is 
it true that the elhptic Fekete polynomials give optimal start equations 
in case n = 1? The experimental data in Figures [3] and H] provides no 
clear support to neither negative nor positive answer. 

The phenomena that are observed in the experiments with one equa- 
tion of degree 4 are amplified when the degree is increased. Figure [5] 
contains the results for degree 10 which show that the random search 
(according to the condition number) gives much better performance 
both for the certified and heuristic homotopy tracking algorithm. (We 
used tighter than default tolerances in the heuristic runs in order to 
get the number of steps closer to that of the certified algorithm. See 
the corresponding Macaulay2 script for the exact settings used.) 





91 


92 


53 


54 


55 


9total 




2.02342 


2.02772 


2.03071 


2.03198 


2.0452 


7.15542 


certified 
heuristic 


11342 
2451.49 


11395.2 
2475.89 


11351 
2466.04 


11393.6 
2441.97 


11258.8 
2437.07 


17737.2 
3376.32 



Figure 5. Search by condition number: n = 1, d = 10. 



4.2.2. Case n = 2. One can show that in order to construct the best 
conditioned total-degree start system we should set r ^ 0.746119 in 
(12. 2p . It is possible to express r in radicals: a computation in Mathe- 
matica determines that r has to be a root of r^(l + 4r^) = 1. 

Figure|6]suggests that better start systems can be found via a random 
search. 





51 


52 


53 


54 


55 


9total 


/^(5) 


1.91848 


1.92438 


1.93129 


1.93511 


1.9386 


2.23607 


#steps 


1259 


1260.72 


1265.15 


1265.33 


1283.8 


1301.15 



Figure 6. Search by condition number: n = 2, di = 
d2 = 2. 

For this case already the search by average complexity becomes quite 
expensive with the current implementation of the certified tracking 
procedure and the computation of the size similar to that in case n = 1 
and d = 4 produces systems much further from optimal than those in 
Figure [61 
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As in the case n = 1, larger examples exhibit larger gap in average 
complexity between the near-optimal systems found with the condition- 
number search and the total-degree start system. In case of two quartic 
equations, n = 2 and di = d2 = 4, the results are shown in Figure [71 





91 


92 


93 


9i 


95 


9total 


Kg) 


4.18566 


4.18963 


4.204 


4.32529 


4.35056 


4.91876 


certified 
heuristic 


17003.6 
8872.14 


17342.5 
8782.57 


17143.8 
8863.36 


17358.6 
8944.19 


17209.6 
8840.86 


20083.3 
9599.85 



Figure 7. Search by condition number: n = 2, di = 
^2 =4. 

We draw the following conclusions: 

• There is a way to find start systems that perform better than 
the total-degree start systems of the form (12. 2p by a random 
search. 

• While bounding the condition number n{g) of a system g can be 
translated into good theoretical bounds (see [3 Theorem 3.7]) 
on the average complexity, the latter does not depend solely on 

• If a system with the optimal fi{g) is explicitly known (as in the 
case n = 1 and li = 4), it is not clear whether it is optimal. 

5. Discussion 

In this article we have demonstrated the possibility of finding start 
systems that lead to better average complexity of finding all roots of 
a polynomial system than the systems of the total-degree family (12. 2p . 
Finding optimal start systems is a hard theoretical problem; however, 
good approximate solutions to this problem can be either constructed, 
e.g., by finding or approximating Fekete polynomials in case n = 1, 
or searched for experimentally. Once near-optimal start systems are 
available, these could be exploited in practice, e.g., in heuristic algo- 
rithms. 

Here are some open questions, the answers to which we would like 
to know: 

(1) If a start system g has the optimal condition number fi{g) , 
does this imply that the average complexity of finding all roots 
(of a target system) with the corresponding linear homotopy is 
optimal? (This is similar to Conjecture [2] that deals with the 
case of just one root.) 
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(2) While the answer to the first question may turn out to be pos- 
itive, one can show that fj^{gi) < ^^{92) does not force the same 
relation for the respective average complexities. However, what 
can be said if the relation between the condition numbers is 
uniform for all roots? I.e., if ^{gi, Zi) < fi{g2, Z2) for all roots 
zi e fi'r^(O) and Z2 G 92^(0), is 

J2 E(Co(/,^?i,zi))< Yl E(Co(/,^?2,;^2))? 

(3) The near-optimal start systems we can find by a random search 
are not optimized for evaluation, which may be a bottleneck in 
the practical computation. Can one construct families of start 
systems that are, on one hand, easy to evaluate and, on the 
other hand, give better average complexity than the currently 
used families of total-degree start systems? 

The capabilities of the random search beyond cases with a small 
number of equations with a modest number of solutions are rather 
limited. Roughly speaking, this is due to the decreasing probability 
of picking a near-optimal system on the sphere S as the dimension 
of S grows. In view of this, a computationally inexpensive answer to 
question (3) above would be particularly valuable. 
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